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B Physics on the Lattice: A, Ai, 7715(7715), A2, 5° — B° mixing, fs and all 
that. 
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We present a short review of our most recent high statistics lattice determinations in the HQET of the following 
important parameters in B physics: the B-meson binding energy, A and the kinetic energy of the b quark in the 
B meson, Ai, which due to the presence of power divergences require a non-perturbative renormalization to be 
defined; the MS running mass of the b quark, mt,(m(,); the B*-B mass splitting, whose value in the HQET is 
determined by the matrix element of the chromo-magnetic operator between B meson states, A2; the B parameter 
of the B°-B° mixing, Bb, and the decay constant of the B meson, /s. All these quantities have been computed 
using a sample of 600 gauge field configurations on a 24 3 x 40 lattice at /3 = 6.0. For A and rfib(mb), we obtain 
our estimates by combining results from three independent lattice simulations at f3 = 6.0, 6.2 and 6.4 on the same 
volume. 



1. The calculation of A and m6(m&). 

The presence of renormalon singularities in the 



divergences. Consider the two-point correlation 
function C (t) for large values of t 



pole mass, m^ ole , implies that the intuitive def- 
inition A = rriB — m^ ole does not correspond to 
a physical quantity. Thus, if A is to be used in 
phcnomcnological applications, a definition free 
of renormalon ambiguities has to be given. One 
can define the B-meson binding energy from the 
experimentally measured value of some physical 
quantity, for example the inclusive B — > XI i>i 
lepton spectrum. It is important to notice that 
such a definition, although correct, depends on 
the order of the perturbation series that deter- 
mines the theoretical prediction for the experi- 
mental quantity. 

Lattice simulations of the HQET provide the 
opportunity to compute A non-perturbatively by 
using the definition of A proposed in ref. jlj which 
is free of renormalon ambiguities. On the lattice, 
renormalons are absent but matrix elements of 
lattice operators may contain power ultraviolet 
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C{t) = J2(o\J B (x,t)4(o,o)\o) 

X 

-> Cieexp(-£U) 

where jjj creates a B-meson from the vacuum, 
is tempting to define A = £. However, the "bare" 
binding energy £ cannot be a physical quantity 
because it diverges linearly as a — > 0, a being the 
lattice spacing |§] . The origin of the linear diver- 
gence is the mixing of the operator hD^h with the 
lower dimensional operator hh. As pointed out in 
J2j , the linear divergence must be subtrated non- 
perturbatively in order not to reintroduce renor- 
malon ambiguities. This divergence can be elim- 
inated from any correlation function by defining 
a renormalized action hDfh of the form M 



hDfh 



hD^h — (1 — exp(— aSm)) hh 
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where 8m is a mass counter-term which can be 
chosen in many different ways. Our preferred def- 
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inition of dm is S 

— 5m = lim — In [ — -pzjzr-r ^ 

{->oo a \ 0(2;, t) 



(3) 



where <5'(x, i) is the heavy-quark propagator in 
the landau gauge which can be computed in a 
numerical simulation of the HQET. In words, (|3|) 
is the condition that the subtracted heavy-quark 
propagator has no exponential fall at large times 
or equivalently that (h(p4 = 0) | hDfh \ h(p 4 = 
0)) = 0. Due to confinement effects, it is not clear 
that @ has a finite limit as t — > 00. Our nu- 
merical results obtained from three independent 
high-statistics simulations at j3 = 6.0, 6.2 and 6.4 
support our assumption on the infrared bahaviour 
of the heavy-quark propagator: no confining ef- 
fects are visible up to distances of about 1.5 fm, 
which is the largest distance which we can reach 
H,|| . An alternative definition of Sm which does 
not require any assumption about the large time 
behaviour of the heavy quark propagator is given 
in refs.|^,Q. However, with this definition A is 
not a parameter of O(Aqcd)- Notice also that 
the divergent part of Sm must be gauge invariant. 
Moreover, the gauge dependence of possible finite 
parts is eliminated from physical quantities up 
to the order of perturbation theory at which the 
calculation is performed. The renormalized bind- 
ing energy A which is of order O(Aqcd), free of 
renormalon singularities and independent of the 
renormalization scale is defined as: A = £ — 5^. 
Our best estimate of A is A = 180±2q MeV ( see 
ref.J| for details). 

We can also compute the MS running mass of 
the b-quark, m(,(?7T(,). By comparing the heavy- 
quark propagator in the continuum with the one 
on the lattice, we obtain H 



m b (m b ) = [M B - A) x C m (a s ,Sm) 
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where C m is the QCD-HQET matching constant 
(see ref.g for details). Taking M B = 5.278 GeV, 
we find m b (m h ) = 4.15 ± 0.05 ± 0.20 GeV §, 
where the first error is statistical and the second 
is an estimation of the error in the truncation 
at 0(a s ) of the perturbative factor C m , derived 
from the bubble-resummation approximation |j] . 
This error can be reduced to about 100 MeV by 



computing perturbative corrections at two-loop 
order. This calculation requires the evaluation 
of the heavy quark propagator in the continuum, 
already done, and on the lattice at 0(a 2 s ). The 
Langevin approach of ref.Q can be used to per- 
form the computation on the lattice. 

2. The calculation of Ai. 

Like A, Ai has no physical meaning since it de- 
pends on the definition that one decides to adopt 
for the kinetic operator. On the lattice, ma- 
trix elements of this operator contain linear and 
quadratic divergences Q . The finite kinetic oper- 
ator, which is free of power divergences, has the 
form 
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hDlh = hD 2 h - — hDfh -%hh 
a a 

where c\ and C2 are subtraction constants which 
we define by imposing the " physical" non-pertur- 
bative renormalization prescription 

(h(p = 0) I hD\h I h(p = 0)) = (6) 

in the landau gauge. This is equivalent to the 
condition pg 2 (t) = c\ + C2 t where 
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with x — (x,t) and y = (y,f). To determine c\ 
and C2, we fit pg 2 (t) for large t to a straight line. 
Our numerical results show that the infrared limit 
of the ratio in eq.(Q) does indeed exist, in spite of 
the confinement effects in the heavy-quark prop- 
agator. C2 can also be obtained directly by elimi- 
nating the sum over t' in eq. (0) and searching for 
a plateau in t and t'. We find that, after sub- 
tracting the contribution of contact terms, both 
methods give indistinguishable results (see ref . 
for details). Notice also that, for the same reason 
as for Sm, C2 must be gauge independent. 

The renormalized kinetic energy Ai is then 
given by: X 1 = Zg 2 (\\ are - c 2 ), where \\ are = 

(B I hD 2 h I B)/(2M B ) and Zg 2 is a gauge in- 
variant matching constant known up to one loop 
order 0] . A convenient way to extract \^ are is to 
study the large time behaviour of the ratio of the 
two-point and the three-point, with an insertion 
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of the kinetic operator, B-meson correlation func- 
tions. The significant reduction of the statistical 
errors with respect to our previous study allows 
for a much better identification of the plateau 
region, which could not be really observed in 
ref.§. We find Xi{B d ) = 0.09 ± 0.14 GeV 2 and 
for the physical parameter Xi(B s ) — \\(Bd) = 
-0.09 ± 0.04 GeV 2 (see ref.g for details). In 
order to make a meaningful comparison of our 
results with other theoretical determinations, a 
perturbative calculation of the relation between 
the different definitions of Ai, and also of A, is 
needed at a sufficient degree of accuracy || . This 
calculation is missing to date. 

3. The calculation of A2 

The chromo-magnetic operator that enters in 
the calculation of A2 is only logarithmically diver- 
gent in the ultraviolet cut-off so that there is no 
need of a non-perturbative subtraction. More- 
over, A2 has a physical meaning because it cor- 
responds directly to the measurable B*-B mass 
splitting up to higher power corrections in 1 / mj 

bare 
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where Z~ is the renormalization constant neces- 
sary to remove the logarithmic divergence present 
in the bare operator. The procedure to extract 
\ h £ re is the same as for \\ are (see ref.0] for de- 
tails). We find A 2 = 0.07 ± 0.01 GeV 2 corre- 
sponding to a mass splitting which is about one 
half of the experimental value. In our opinion, 
there are three possible explanations of the dis- 
crepancy, which is common to all lattice results. 
The first is the fact that, at one loop order, Z g .Q is 
very large indicating that higher-order terms may 
modify this factor significantly. A second source 
of systematic error is quenching. The final pos- 
sibility is that the discrepancy is due to the fact 
that the mass splitting receives significant contri- 
butions from higher terms in the HQET. Further 
investigation is needed on this subject. 

4. The calculation of Bb and Jb 

We have also computed the parameter f B Bs, 
which enters the theoretical prediction of the B°- 
B° mixing parameters Xd and x s . We measure Jb 



and Bb from the ratios of three- and two-point 
functions with the AB = 2 operator fixed at the 
origin. Due to the high statistics we have accu- 
mulated, we can study this ratios at large time 
distances at which the ground state has been un- 
ambiguosly isolated. We find a very good signal 
for all the correlation functions and ratios. How- 
ever, the main uncertainty in our results comes 
from the evaluation of the renormalization factors 
relating the effective operators on the lattice to 
those in the continuum. For Bb, we use the new 
calculation of the NLO Wilson coefficients of the 
0(a s ) full theory-HQET matching, which takes 
into account previously missed contributions 
In our opinion, a non-perturbative computation 
of these factors is crucial to reduce the rather big 
systematic errors. This work is in progress now 
|J. Our best estimates are @: B Bd = I.21±0.06 
(RGI B parameter), B B jB Bd = 1 .011 ± 0.008 
and I 2 B s B B JpB d BB d = 08 ±0.07, where we use 
the boosted perturbation theory to evaluate the 
renormalization constants. 
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